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1. Introduction 

Our use of lattice field theory to calculate the structure of the nucleon from first principles has 
two complementary objectives. One goal is to achieve quantitative agreement with experimental 
observables such as nucleon form factors and parton distributions, to confirm the quantitative pre- 
cision of our solution of QCD and to establish the credibility to make predictions and guide future 
experiments. However, merely producing a black box that only reproduces experiment would be 
unfulfilling. Hence, our second goal is to obtain insight into how QCD works, revealing, for ex- 
ample, the origin of the nucleon spin, physical mechanisms, such as instantons, responsible for 
essential features of hadron structure, and the dependence on parameters of QCD like N& Nf, m q , 
and the gauge group. A crucial issue in making quantitative contact with experiment is calculating 
sufficiently far into the chiral regime that reliable chiral extrapolations to the physical pion mass are 
possible. Hence, in this work we have utilized the extensive set of configurations with dynamical 
improved staggered quarks generated by the MILC collaboration ^ to explore nucleon structure 
in the chiral regime. 

2. Mixed Action Lattice Calculation 

As explained in Ref. [0], we utilize a hybrid action combining domain wall valence fermions 
with improved staggered sea quarks. Improved staggered sea quarks offer the advantage that due to 
the relative economy of the algorithm, lattices with large volumes, small pion masses, and several 
lattice spacings are publicly available from the MILC collaboration. Although the fourth root of 
the fermion determinant remains controversial, current evidence suggests it is manageable [Q, H]. 
Renormalization group arguments indicate that the coefficient of the nonlocal term approaches 
zero in the continuum limit partially quenched staggered chiral perturbation theory accounts 
well for the artificial properties at finite lattice spacing and the action has the advantage of 
being improved to &(a 2 ). Domain wall valence quarks offer equally compelling advantages to 
justify investing resources in calculating hadron observables on staggered configurations that are 
roughly comparable to the resources required to generate the configurations themselves. Domain 
wall fermions prevent mixing of quark observables by chiral symmetry, are accurate to 0(a 2 ), and 
possess a conserved five dimensional axial current that facilitates calculation of renormalization 
factors. In addition, hybrid action (often referred to as mixed action) chiral perturbation theory 
results are available for many observables, and by virtue of an exact lattice chiral symmetry, one 
loop results have the simple chiral behavior observed in the continuum. The parameters of the 
configurations used in this work are shown in Table |[lj], and details can be found in Refs. [§, 0]. 

3. Moments of Parton Distributions 

Parton distributions measure forward matrix elements of the gauge invariant light cone opera- 
tors 

ff r (x) = / —e iXx q(-Xn I2)Y 3»e- ig ^ 2danA{an) q(Xn/2) , (3.1) 
J 4n 

where x is a momentum fraction, n is a light cone vector and r = tf. or Y = rf,y^ . Using the operator 
product expansion, the operators in Eq. [ |3.1| ] yield towers of symmetrized, traceless local operators 
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Table 1: Lattice parameters used in this work. 



that can be evaluated on a Euclidean lattice 

= qy {tll [r5}iD^---iD^q, (3.2) 

where [75] denotes the possible inclusion of 75, the curly brackets represent symmetrization over 
the indices pLj and subtraction of traces, and D= 1/2(D — D). A related operator for transversity 
distributions is 

= qy 5(J ^HD^---iD^q. (3.3) 

Using the notation and normalization of Ref. the forward matrix elements (P,S\^^ 1 '" lln+1 ^\P,S) 
yield moments of the unpolarized quark distribution: 

(x") q = f l dxx n [q(x) + (-l) n+1 q(x)], (3.4) 
JO 

the forward matrix elements (P,S\ffy^ l '"^" +1 ^\P,S) yield moments of the helicity distribution: 

(x") Aq = [ dxx"[Aq(x) + (-l) n Aq{ X )], (3.5) 
Jo 

and the forward matrix elements (P, S\ iff^ 1 " ^ \P, S) yield moments of the transversity distribu- 
tion: 

(x n ) Sq = J dxx n [dq(x) + (-l)' I+1 5^(x)]. (3.6) 

In this work, we calculate only connected diagrams, and hence concentrate as much as possible on 
isovector quantities. 

All quark bilinear operators in Eqs. [3^] and [3_3] are renormalized as follows The axial 
current is renormalized exactly using the conserved five dimensional axial current. By virtue of the 
suppression of loop integrals by HYP smearing, the ratio of the one-loop perturbative renormal- 
ization factor for a general bilinear operator to the renormalization factor for the axial current is 
within a few percent of unity, suggesting adequate convergence at one-loop level. Hence the com- 
plete renormalization factor is written as the exact axial current renormalization factor times the 
ratio of the perturbative renormalization factor for the desired operator divided by the perturbative 
renormalization factor for the axial current. 
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3.1 Chiral Perturbation Theory 

Ideally, we would like to perform high statistics calculations at pion masses below 350 MeV 
and extrapolate them in pion mass and volume using a chiral perturbation theory expansion of 
sufficiently high order to provide a quantitatively controlled approximation. In practice, our most 
convincing chiral extrapolation has been for gA using the finite volume results including A interme- 
diate states of Ref. JTOj], where the fit involving 6 low energy parameters yielded an excellent fit up 
to the order of a 700 MeV pion mass and agreed with experiment with 6.8% errors [O]. Similar 



extrapolations of gA have been performed by other groups [12]. This success for g A is particularly 
relevant to the subsequent discussion of nucleon spin, because it involves the same operators (l) q 
as AZ. We note that because the nucleon and A should be included together at large N c [13] and 
indeed show large cancellations in the axial charge, we prefer to include the A as an explicit degree 
of freedom in the analysis. 

An unresolved puzzle in calculating moments of structure functions is the relatively flat behav- 
ior of the momentum fraction (jc) at a constant value substantially higher than experiment [14, 0]. 
Hence, it is particularly interesting to ask whether a chiral perturbation theory fit determined with- 
out knowledge of the experimental result is in fact statistically consistent with experiment. Since 
there is presently insufficient data to perform a full analysis including the A, here we present a sim- 
ple self-consistent improved one-loop analysis using only nucleon degrees of freedom that appears 
to work very well in our regime. 

The details will be presented in a future publication [|lg], but the basic idea is as follows. We 
begin with the one loop expression at scale /I [|i7,|Iq]. 

(3«lo + l) '~ 2 " ' 



{X n )u-d 



in which we explicitly note that gA$ and f na are g A and f n in the chiral limit. We are free to choose 
the scale \i to be f n . Additionally we replace gA.o and f nfi with their values at the given pion mass 
gA,m K and fnj7i K , so that the result may be rewritten as 



(x n ) u _ d =a n 1 ~ T\' ' "** In l-p- +K^. (3.8) 




If we view this as an expansion in the ratio r = , . , ; % — , one can show that shifting to an expansion 
around gA and f n defined at another mass only introduces changes of &{r 2 ). Hence, to leading 
order, we may write an expression in which we use the values gA,iat, f%,lat, and m n ^ at calculated 
on the lattice at specific values of the quark mass. Then, the expressions for the moments of the 
unpolarized, helicity, and transversity distributions are the following: 



{x )u-d = a n ( 1 jf-^ - 7 — In | -pj— ) ) +b n —^- (3.9) 

J JT.lat 




"4,lat 



Km 



(^)Au-M = Afl» ( 1 - 7;^ 2 ,2 h I I I +Mn-^ (3.10) 

/^,lat 



"4,lat 



^)su- 5d = Sa n \ 1- Z'TL, J? 1° I 1 \+8b n -^. (3.11) 
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Figure 1: Zeroth moment of helicity distribution. Figure 2: Zeroth moment of transversity distribution. 




Figure 3: First moment of unpolarized distribution. Figure 4: First moment of helicity distribution. 



These results allow a least-squares two-parameter fit to the lattice data for moments and provides 
an extrapolation to the physical pion mass with a corresponding error band. Note that the series 
is substantially rearranged, by virtue that the calculated values of gA, fn, and m K are used at each 
value of the bare quark mass. Although we cannot prove that this self-consistent improved one-loop 
result should be accurate throughout the range of our data, to the extent to which it is successful, we 
believe its success arises from this self-consistent rearrangement. Additionally, the use of physical 
rather than chiral limit values for f K was first tried in [19] and has since been studied in chiral 
perturbation [20, 21] and applied to a variety of lattice calculations [22, 24, 25, 59] . 



3.2 Lattice Results 

Here we show the results of this one-loop analysis. Figure [Jl]] shows the result for gA, 
which is nearly as good as the complete analysis of Ref. [|j]], and yields a comparable extrap- 
olation and error bar. Reassured by this result, we show analogous results for (l)s u -Sdr ( x )u-d, 
(x)Au-Ad,(x)8u-8d, and ( x2 )su-Sd in Figs. [g]-[@|. 

Note that in every case for which there is experimental data, this analysis, which in no way 
includes the experimental result in the fit, yields an extrapolation consistent with experiment. The 
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Figure 5: First moment of transversity distribution. Figure 6: Second moment of helicity distribution. 



results are collected together in Fig. [Q], where because experimental results are not available for 
all cases, we have normalized all results to the corresponding lattice result. 
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Figure 7: The six moments considered in this work. Lattice results are shown in blue and experimental 
measurements in red, and each is normalized to the corresponding lattice result. 



4. Form Factors 

Form factors are interesting physically because at low momentum transfer they characterize 
the spatial size of charge and current distributions and at high momentum transfer they measure 
the ability of the nucleon to absorb a large momentum and distribute it to all the constituents such 
that the system remains in its ground state. Although for a relativistic system, the slope of the 
form factor is not precisely related to the rms radius, we will adhere to the common usage and 
refer to the slope as the rms radius. (We note in passing that the slope of F\ is in fact related to the 
transverse rms radius in the infinite momentum frame.) Our primary focus here is on the qualitative 
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approach to experiment as the quark mass is decreased, and to avoid uncalculated contributions of 
disconnected diagrams, we consider isovector form factors. 

The nucleon vector current form factors F\ and F2 are defined by 

{pWY^W) = u( P )[F l (q 2 )r+Hq 2 ) j ^ L W). (4-1) 

2m 

Figure [||] shows the lattice data and dipole fits for F\ at five pion masses, and one observes that 
the lattice results systematically approach the experimental curve as the pion mass decreases to 
359 MeV. One can quantitatively observe how the rms radius (r 2 )" - '' defined from the slope ap- 
proaches the experimental value by fitting it with the simple chiral extrapolation formula [27] 

{r } =a - log , (4.2) 

(47TA) 2 \m 2 K +A 2 J 

where A is a phenomenological cutoff or equivalently, coefficient of an m\ term in the expansion. 
Note that, in contrast to most chiral extrapolations which contain finite terms of the form m\ logm 2 , 
the isovector radius diverges like log/n 2 , rendering the variation of the radius quite substantial near 
the physical pion mass. Figure ^ shows the results of fitting the lattice data with Eq. [4- .2], and 
one notes that without having been constrained to do so, the chiral extrapolation is consistent with 
experiment. 

The form factor F2 is of particular interest following the observation at JLab p9| ] that measure- 
ment of spin transfer yields a form factor that decreases more slowly with momentum transfer than 
the traditional Rosenbluth separation, which is now believed to suffer from substantial contami- 
nation from two photon exchange contributions. Figure [ |l0[ ] shows that the lattice results indeed 
approach the experimental ratio FijF\ from spin transfer as the pion mass decreases. 

The nucleon axial vector current form factors Ga and Gp are defined by 

(pWfKY\p') = u(p)[G A (q 2 )fY 5 + ^y5G P {q 2 ) + a^y 5 q v GM{q 2 )W). (4.3) 



Figure [11] shows lattice calculations of the isovector ratio Ga/F\ compared with experimental 
form factors extracted from pion photoproduction and neutrino scattering. Again, to within the 
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Figure 10: Isovector form factor ratio F%lF\ at three Figure 11: Isovector form factor ratio Ga/F\ at four 
masses compared with experiment [E8|. masses compared with experiment. 



0.2 0.4 



¥ Expt: electroproduction 
x Expt: muon capture 

• m ]t = 775MeV 

♦ m* = 605MeV 
m % = 498 MeV 

■ m K = 359 MeV (2.5 fm 3 ) 

° m K = 359 MeV (3.5 fm 3 ) 
► m =301 MeV 



0.6 0.8 
Q 2 (GeV 2 ) 



1.2 



Figure 12: Isovector form factor ratio Gp/G^ at six Figure 13: 
masses compared with experiment. malized as 



Isovector pion pole contribution, nor- 

q 2 +m 2 K )Gp/4m 1 G A , at m n =350 MeV 



roughly 10% discrepancy between the two experiments, the lattice results are qualitatively consis- 
tent with experiment. In the soft pion limit, Gp is dominated by the pion pole: 



2 4m 2 G A (q 2 ) 
Gpiq } ~ qi+ml ■ (4 " 4) 



Figure [12] shows lattice results for the the ratio Gp/Ga compared with experiment. To demon- 
strate the degree to which Eq. [[Oj] is satisfied, Fig. [ ]T3[ ] shows the ratio (q 2 + m 2 z )Gp/4m 2 GA for 
lattice data at pion mass 350 MeV, which is consistent with unity. 

5. Generalized Parton Distributions 



Off diagonal matrix elements of the tower of operators in Eq. [ p.2| ] yield the generalized form 
factors 

(P'\^\P) = ((^))A 1Q (t) + ^-{((J^ a ))A a B w (t), 
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Figure 14: Generalized form factors A 20, #20 and C20 for u — d and u + d at four masses. 

(P^iM2}|P) = (( 7 ^}))A 20 (0 + «<T' fa > a ))A fi2o(0 + -A^>A^C 20 (0 , 

zm m 

(P>\0{^Wi}\P) = p{MipM2(( 7 M3})) A3() (;) + J-p{Mip/fa(( <y M3}«))A a B3o(0 

2m 

+ A^A^^(Y^)A 32 (t) + ^A^A^((o^ a ))A a B 32 (t), (5.1) 

where we use the short-hand notation ((r)) = U(P',A')rU(P,A) for matrix elements of Dirac 
spinors U and where A = P' — P and t = A 2 . 

Of particular interest in this work is the relationship between the generalized form factors and 
the origin of the nucleon spin. The contribution of the spin of the up and down quarks to the total 
spin of the nucleon is given by the zeroth moment of the spin dependent structure function (1)^ as 

^AZ=^{\) Au+Ad . (5.2) 
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Figure 15: Nucleon spin decomposition. Squares Figure 16: Nucleon spin decomposition by flavor, 
denote A£ M+d /2, the star indicates the experimental Squares denote AE"/2, diamonds denote AE d /2, tri- 
quark spin contribution, and diamonds denote L u+d . an § les denote L "' and circles denote L<i - 



Note both that our previous calculation of gA confirms our ability to calculate the connected con- 
tributions to (I)a^ accurately on the lattice and that A£ requires the calculation of disconnected as 
well as connected contributions. Throughout this section, we will only discuss the results of con- 
nected diagrams, so all results will eventually need to be corrected for the effect of disconnected 
diagrams as well. 

The total contribution to the nucleon spin from both the spin and orbital angular momentum 



of quarks is given by the Ji sum rule [ 3C ] 



1 



^=2( A 20 d (0)+*20 +d (0) 



(5.3) 



so that the contribution of the orbital angular momentum is given by L q = J q — jA£. Earlier calcu- 
lations in the heavy quark regime showed that for 700 MeV pions, roughly 68% of the spin of the 
nucleon arises from the spin of quarks, 0% arises from orbital angular momentum, and hence the 
remaining 32% must come from gluons [31, 32, 33]. Figure [IP]] shows the recent lattice data for 



A20, #20 and C20 for lighter pion masses, and the results for AT and L q are shown in Fig. [ 15]. The 
full decomposition showing the contribution of spin and orbital angular momentum from up and 
down quarks is shown in Fig. Jlrj]. Here, one observes that the angular momentum contributions 
of both up and down quark are separately substantial, and it is only the sum that is extremely small. 



In previous work in the heavy quark regime []34|], we have emphasized the dramatic differences 
in the slopes of A 10, A20, and A 30 and the fact that this reflects a sharp decrease in the transverse 
size of the nucleon as x approaches 1 . To show that this behavior also arises in the chiral regime, 
Fig. [ |l7| ] shows that the slope of the form factor ratio A30/A10 differs substantially from unity as 
the pion mass is decreased and also approaches the result given by a phenomenological parameter- 



ization of the generalized parton distributions [35]. 



6. Conclusions 

In summary, the hybrid combination of valence domain wall quarks on an improved staggered 
sea has enabled us to begin to enter the era of quantitative solution of full lattice QCD in the chiral 
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Figure 17: Comparison of the ratio A30/A 10 for u — d and u + d at four masses with a phenomenological fit 
to generalized parton distributions. 



regime. The axial charge, g&, represents a successful, "gold-plated" test, which argues well for the 
prospects of quantitative control of a range of important nucleon observables. The chiral extrapo- 
lation of moments of quark distributions using our self-consistently improved one-loop analysis is 
encouraging, but of course we would still like to directly calculate the turn over in the approach to 
the chiral regime. Similarly, the general agreement between nucleon form factors of the vector and 
axial currents is highly encouraging. Generalized form factors are also being calculated well into 
the chiral regime, and give promise for understanding the origin of the nucleon spin and the trans- 
verse structure of the nucleon. In the long term, since lattice calculations determine moments of 
generalized parton distributions and experiments measure convolutions of generalized parton dis- 
tributions, there is an excellent opportunity for synergy between experiment and theory in jointly 
determining generalized parton distributions. The final analysis of all the data shown in this work 
is currently being completed, and full results will be published in the near future. 
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This work also indicates obvious challenges for the future. Clearly the calculations must be 
extended to lower quark masses and finer lattices and analyzed with partially quenched hybrid chi- 
ral perturbation theory, and disconnected diagrams must be calculated. In addition, it is important 
to develop new techniques to explore form factors at high momentum transfer, gluon observables, 
and transition form factors for unstable states. 
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